Anderson localization in a periodic photonic lattice 
with a disordered boundary 
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We investigate experimentally the light evolution inside a two-dimensional finite periodic array of weakly- 
coupled optical waveguides with a disordered boundary. For a completely localized initial condition away from 
the surface, we find that the disordered boundary induces an asymptotic localization in the bulk, centered 
around the initial position of the input beam. © 2012 Optical Society of America 
OCIS codes: 130.2790, 240.6700, 290.5825 



The original concept of Anderson localization assumes 
a periodic structure where disorder is introduced via a 
random change of the local properties at each site of the 
lattice, leading to wave localization due to interference 
between multiple scattering paths [1]. Although it was 
first described in the context of condensed matter [1,2], 
there are now examples in many other fields: Acoustics 
[3], microwaves [4], Bose-Einstein condensates [5] and 
optics [6-10], to name a few. Most importantly, in the 
optical domain it was proven in experiment and theory 
that a random displacement of the lattice sites yields the 
same localization results as changing the local properties 
of the individual sites [11, 12]. 

In rough or corrugated channels, the phenomenon of 
Anderson localization is connected to the transmission of 
electrons or optical pulses. In these systems, the chan- 
nel surface is disordered along the propagation direc- 
tion. Theoretical studies showed transitions between dif- 
fusive and localized regimes [13-15]. Transport behavior 
was also considered in [16-18], where multiple scatte- 
ring from longitudinal surface roughness caused localiza- 
tion of waves, although regimes of coexistence between 
ballistic, diffusive and localized transport, depending on 
the symmetry of the corrugation profile, have been pre- 
dicted [19]. 

Recently, a different kind of "corrugated waveguide" 
was considered theoretically [20]. Here the disorder is 
only in the transverse direction, and does not change 
along the direction of propagation. Also, the optical 
medium in the bulk, away from the corrugated surface, 
possesses a periodic index of refraction along the two 
transversal directions, thus forming an array of weakly 
coupled waveguides [21]. The disorder is imposed on the 
boundary by random displacement. Despite the weak 
disorder, a light beam still tends to localize in the center 
of such a system far from the boundary. In our work, we 
experimentally prove this prediction. 

Let us consider a 2D rectangular N x M waveguide 




Fig. 1. Scheme of a finite two-dimensional coupled array 
of elliptical waveguides with disordered boundary. 

array (Fig.l). In the coupled-mode approach, the elec- 
tric field E{r, z) propagating along the waveguides can 
be written as a superposition of the waveguide modes, 
u{r,z) = X]n'^n(^)0(r — n), where r = (x,^), is the 
amplitude of the single modes 0n centered around site 
n = (p, q) . The evolution equations for the mode ampli- 
tudes Un are 

Yl Yl ^hjUpj = -i^, (1) 

j=P±i j=q^'i- 

where n = (p, q) denotes the position of the guide cen- 
ter, z is the longitudinal distance, and the Cy^h are the 
coupling coefficients between nearest-neighbors guides. 
They decay exponentially with the mutual distance be- 
tween the guides. To keep our approach general, we as- 
sume anisotropic coupling in the horizontal and vertical 
direction. At the boundary of the array, randomness is 
introduced by displacing the guides from their (ordered) 
positions, along the boundary surface. This creates ran- 
dom couplings among the boundary guides, for the hori- 
zontal coupling Ch Che~^^^ as well as for the vertical 
coupling Cy Cye~^'"^. The quantity A is a random 
number in [—1, 1] and the randomness strengths ujh and 
ujy along the horizontal and vertical boundaries respec- 
tively are different due to the ellipticity of the guides. 
The bulk, that represents the ordered part of the array, 
is connected to the disordered boundary layer with cou- 
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pling values computed from the mean pythagorean dis- 
tance. These conditions are close to the possibilities of 
experimental realization [12]. The entire system could be 
viewed as a single "corrugated" photonic crystal waveg- 
uide, with fixed transversal corrugation. 

First, in order to illustrate the localization mecha- 
nisms in the system under consideration, we analyze a 
waveguide array with 13 x 5 waveguides and disordered 
boundary and perform a direct numerical integration of 
Eq. (1) with a single-site excitation at the array center 
{n Carrie) = (7,3). The parameters chosen for the sim- 
ulation of the disordered case (left column of Fig. 2) 
are: Ch = 2, = 1.7, Uh = 11/28, ujy = 4/17 [22]. 
For comparison, we compute the evolution for the or- 
dered case and its pure discrete diffraction (right col- 
umn of Fig. 2). For small propagation distances {z 5)^ 
both systems show discrete diffraction with maxima at 
the outer lobes. The amplitude of the initially excited 
site decreases steadily with z. However, after a transi- 
tion distance of Zt — 10, that is necessary for a com- 
plete backscattering cycle, the infiuence of the boundary 
layer becomes apparent. In the case of the disordered 
boundary, for z > 10, the amplitude of the center site 
is non- vanishing for all propagation lengths. It continues 
to oscillate with z, but finally saturates to a nearly sta- 
tionary mode, showing the persistence and stability of 
localization due to the disordered boundary. The main 
difference to the case of a completely disordered system 
(bulk + surface) is that the initial oscillations due to 
discrete diffraction are stronger. 

For our experiments, we fabricated various waveguide 
arrays in polished bulk fused silica glass, using the laser 
direct- writing technology [23]. Each guide has dimen- 
sions of 4 X 12 /im^ and exhibits a refractive index change 
of ~ 5 X 10~^ [22]. We prepared 60 disordered-boundary 
waveguide arrays with a length of 101 mm, with either 
ArxM = 5x5or with x M = 13 x 5 waveguides each. 
The inter-guide separation in the ordered bulk, which 
is also the mean separation at the (disordered) bound- 
ary, was of 14 fim and 17 /im, along the horizontal and 
vertical directions, respectively. Disorder in the waveg- 
uide spacing was induced by varying the inter-guide dis- 
tance: dh = 5.5A /im and dy = 17 ± 4A /im, with 
A G [—1, 1] randomly equally distributed. 

In each array the individual central waveguide was 
excited using a Ti: Sapphire laser system at low input 
power to ensure linear propagation. At the end facet, 
the intensity patterns were recorded with a CCD cam- 
era. For the observation of localization in the 5 x 5 ar- 
rays, a wavelength of 800 nm was used, corresponding 
to coupling coefficients [22] of Ch{Cy) 2.0(1.7) cm~^. 
To enhance localization in the arrays of 13 x 5 waveg- 
uides, the wavelength was increased to 840 nm, leading 
to higher coupling coefficients and therefore decreased 
coupling lengths. This allows the observation of localiza- 
tion effects at effectively shorter evolution distances. 

The experimental results, shown in Figs. 3 and 4, 
demonstrate a clear localization tendency around the po- 
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Fig. 2. Comparison between simulations of an ordered 
13 X 5 array (right column), and the disordered boundary 
waveguide array) of eq. (1) (left column) with Cq, = 2, 
C(3 = 1.7, uJa = 11/28, cj/j = 4/17, averaged over 100 
realizations. 



sition of the input beam, for both, 5x5 (Fig. 3) and 13x5 
(Fig. 4) arrays. Figure 3(a) shows the output intensity 
at the end facet, averaged over 30 different boundary- 
disorder realizations. The exponential localization be- 
havior of the light intensities in the horizontal direction 
is shown and clearly seen in Fig. 3(b). Since the global 
coupling in the vertical direction is weaker, the decay 
here is remarkably stronger (Fig. 3(c)). 

A microscopic image of a single realization of the 13x5 
disordered boundary array is shown in Fig. 4(a), where 
the ordered center as well as the boundary layer with 
disordered spacing can be seen. For comparison we also 
show a microscopic image of the completely ordered 
13 X 5 waveguide array. In Fig. 4(c), the mean output 
intensity of 30 realizations is shown, where localization 
on the central input site is clearly observed. When com- 
paring with the intensity distribution in an ordered array 
(Fig. 4(d)) we find that the propagation distance is suf- 
ficiently long for complete delocalization in the ordered 
case. Thus, the observed localization effects can be at- 
tributed to the disordered boundary layer. The stronger 
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Fig. 3. Experimental results for the 5x5 arrays: (a) 
mean output of 30 realizations, A = 800 nm, (b), (c) 
logarithmic plots of \un,mj^ \unc,m\ respectively. 



coupling obtained in the 13 x 5 array due to the longer 
wavelength leads to better exponential localization in 
both coupling directions, shown in Figs. 4(e) and 4(f). 

In summary, we have investigated experimentally the 
influence of a disordered boundary in a finite 2D coupled 
waveguide array and found asymptotic partial localiza- 
tion of the wave packet in the center of the bulk region 
far away from the boundary. We conclude that the pres- 
ence of a disordered boundary can give rise to Anderson 
localization in regions away from the boundary. We con- 
jecture that this result could be extrapolated to larger 
finite lattices, if one allows for a sufficiently long propa- 
gation distance. 
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